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We present a first principles account of Bloch Walls on the basis of relativistic Spin-Density 
Functional Theory in the Local Density Approximation (LDA). We argue that calculations based 
on this description will provide useful and novel informations about the physics of this important 
O \ element of micro-magnetics. To illustrate the points we wish to make we implement the proposed 

■ calculations for Fe and determine the Bloch Wall thickness Ibw. Moreover we study features of the 

electronic structure which arise due to the presence of the Bloch Wall. 
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I. INTRODUCTION 



Domains of different magnetic orientations spontaneously nucleate in all magnetic materials in order to minimize 
the magnetostatic energy (dipolar interaction) of the system. These domains are separated by domain or Bloch Walls 
where the magnetization is changing over, rapidly, from one orientation to an other. In this paper we shall develop 
a theoretical account of these transition regions in terms of mobile electrons and their spins, that is to say from first 
^ ' principles. 

The understanding and control of the domain structure, magnetic morphology, is central to the effective use of 
^ ■ magnetic materials in many technologies O. Moreover, the nature, shape, and motion of Bloch Walls have always 
<^ ' attracted fundamental scientific interest (2[^]. Nevertheless, all the theoretical discussions so far, of the entire subject, 
QQ have been strictly phenomenological. In the present paper we wish to lay the foundation of a first-principles approach 
to the problem. 

■ To start with a tractable problem, henceforth, we shall focus on that of a single Bloch Wall. Note that whilst its 
width, ^BW ~ 10 — 1000 nm, much larger than the lattice spacing a, it is smallest of all the other length scales such 
as the sizes of the domains, 1 — 1000 /im, magnetic texture, > 0.1 mm, etc. in the problem. Thus, by deploying our 
first-principles methodology, designed to treat variations in properties on the length-scale a, to Bloch Wall problems, 
we are making the natural choice of crossing the smallest scale-gap first. 
' The established theory of Bloch Walls is a topic in micro-magnetics and it is wholly phenomenological. It consists 
of the equations of Landau and Lifshitz [Q as generalized by Brown For clarity and easy future reference, we 
^ . summarize here the part of this theoretical framework relevant to our present concern. In short, the magnetization 
] density is described by the vector field M{f,t), which evolves in time according to the Landau-Lifshitz equation, 

o: dM^2^.^mM] 

dt h 5M ^ ' 



where /xb is the Bohr magneton, Ti is the Planck constant and is a generalized free-energy functional of M{r, t). 

For our present purposes it will be sufficient to take the .?^[M] to be given by the usual Ginzburg-Landau expansion, 

T = T[M] « j dr\ |VAf.|2 + ^^]M^MJ + ... I (2) 

where i and j refer to Cartesian axes x,y and z, and the coefficients A, Bij are material specific parameters which 
depend on such thermodynamic variables as temperature, T, and pressure, p. Evidently, the equilibrium magnetiza- 
tion, Meq, which minimizes J-"[M] is a stationary solution of Eq. ([|). The simplest Bloch Wall like solution of these 
equations follow from parameterizing the Cartesian components Mi in terms of spherical coordinates 9 and (j) (see 
Fig. 1.) 
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Mx = Mo cos (j)sm9, My = Mq sin (j)sm9 , = Mq cos 9 (3) 

where Mq is a constant equal to the saturation magnetization, and stopping the expansion in Eq. (||) at the 4*'' order. 
For a cubic system {Bij — B 5ij) this means 

T{e, <P)= J {aiV9{r) f + sm^{9){V<j){f)f + ^sin^ ,/)(r) cos^ dr, (4) 

where the new coefficients are defined as a = AMq and /3 = BMq, respectively. Taking 9 — ^ and (p a function of 
the z coordinate only, it can be readily shown |9| that, subject to the boundary conditions (j){~oo) ~ and ^(oo) = ^, 
the function 



4>{z) = arctane 




minimizes ^{9, (f>) in Eq. (^. In what follows we shall study this simple example of a 90° Bloch Wall from the point 
of view of mobile electron whose spin gives rise to such magnetic configuration. 

Given the success and the general satisfaction with the established phenomcnological theory outlined above one 
may, at this stage, protest that such microscopic considerations are not needed. To answer this objection we shall 
now pause briefly to argue the case to the contrary. In the interest of economy we present our arguments as a list of 
brief statements: 

a. ) Firstly, the phenomcnological theory assumes that M (z) rotates as z goes from — oo to oo without a change in 

magnitude. A first-principles theory should tell us if and when this is true. 

b. ) Because, as yet, experiments can not tell us otherwise the phenomcnological theory accepts the form of 

(j){z) in Eq. (|5|) as immutable. By contrast, the microscopic theory could discover situations (materials) where 
this is not so. The same thing can be said about deviations of 9 from ^ as a function of z or, indeed, the xy 
dependence of both 9 and <j). 

c. ) Of course, first-principles calculations not only yield M{r) but also the changes in the electronic structure due 

to the deviation of M (r) from its saturation value, namely the presence of a Bloch Wall. In other words, they 
describe the electronic structure which supports, consistent with, a Bloch Wall. This novel information becomes 
available for the first time in the calculations we shall report. 

d. ) The change in the electronic states due to the presence of a Bloch Wall is, or can be viewed as a description of 

the electron-Bloch Wall interaction. Clearly, a study of this interaction will make an important contribution to 
understanding the widely observed magneto-transport phenomena associated with electrons scattering off Bloch 
Walls 

e. ) Inevitably, a magnetic inhomogeneity will distort the underlying lattice and such distortions give rise to Bloch 

Wall-Bloch Wall interactions and pinning of Bloch Walls by lattice defects like impurities vacancies or grain 
boundaries . Although we do not deal with the appropriate generalization of our first-principles approach, we 
wish to stress that it can be readily adopted to address these issues and this possibility is an important part of 
the motivation for perusing it. 

Finally, we comment on our choice of BCC iron as the host to the Bloch Wall in our calculations. As is well 
known the saturation magnetic moment per atom /Xg at T = and the effective moment deduced from the Curie 
constants, measured at T > Tc, for metallic ferromagnets are usually not the same. In fact, the deviation of fis/fJ-cS 
from 1 can be taken as a good measure of how independent the local magnetic moments are from their relative 
orientation. Since BCC iron is a famously good moment system, namely fig — fJ^cS, we expect the first-principles 
calculations to map optimally onto the phenomcnological theory summarized by Eqs. (^) and This is indeed 
what wc find and we take this fact as a validation of our conceptual framework as well as our numerical procedures 
for implementing it. The prize to pay for having done this 'easy' case is that the importance of all the novel features 
listed above (a.-e.) are minimized. Thus, if we disregard these interesting but small effects, our calculations reduce 
to a new, if complicated, way of computing the material dependent parameters, e.g. the spin-wave stiffness constant 
a and the anisotropy parameter /3. As these are also available from other type of first-principles calculations [p"5|-p^, 
the comparison of our results with those of others in the field is a useful exercise in its own right. 

In the next section we shall describe the first-principles theoretical framework, based on a relativistic, spin-polarized 
density functional theory, for our calculations. This is followed by a section describing the computational procedures. 
We shall present and discuss our results in Sec. 4., whereas in Sec. 5. we shall evaluate the progress we have been 
able to make and the prospects of the microscopic approach we have advocated. 
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II. A FIRST-PRINCIPLES THEORY OF BLOCK WALLS 



A fully relativistic Density Functional Theory (DFT) fl^ill^l , which includes the dipolar interaction between elec- 
trons, would, presumably, yield an inhomogeneous ground state with the domain structure determined by the size and 
shape of the sample. Clearly, the corresponding calculations are out of the question and we shall follow the logic of the 
phenomenological theory. Namely, we neglect the dipolar interaction and study a single Bloch Wall engendered by the 
constraint that the magnetization M(r) is oriented along two different easy axes at z — > ±00. This problem is readily 
encompassed by the spin-polarized relativistic Density Functional Theory in the Local Density Approximation (LDA) 



as usually applied in solid-state physics [[l9y20[ . There are only two features of the way we shall proceed which deserve 
further general comment. Firstly, we note that we shall not be looking for a ground state, as usual, but a lowest 
energy state consistent with a constraint which prescribes a symmetry different from the ground state. Fortunately, as 
it is well known, DFT covers this eventuality Secondly, we shall not attempt to address the vexing conceptual 

difficulties that arise in connection with imposing constraints while solving self-consistently the Kohn-Sham-Dirac 
equation of the theory but in this preliminary exploration of the subject implement our strategy in the simplest 
possible way as outlined below. 

There is no reason to doubt that the established phenomenological theory in the introduction gets the essential 
physics of the Bloch Wall formation right. From the point of view of formulating a first-principles version of this 
theory it may be summarized as follows: the changing orientation of M{z), namely the variation of 4>{z) in Eq. (|^), 
across the Bloch Wall implies an exchange energy cost which is lower the slower the variation, whereas the anisotropy 
energy favors a rapid change from one easy direction, (f> — 0, to the other, = ^, and the width of the transition 

region of 4){z) in Eq. (|^), ^bw — '^\J^-^ determined by the balance of these two tendencies. Indeed, we can bypass 

the solution of the Euler-Lagrange equation, which in the case of minimizing the free-energy J^{6, (p) in Eq. (H) is the 
famous Sine-Gordon equation whose solution is in Eq. (||), and take 0(z) to be a simple function, (/)o(^), which 
goes from to ^ in a distance of L and minimize the corresponding free-energy J'{L) with respect to the width L. 
As can be readily shown, by substituting (j)o{ j;) into Eq. (|^), for ^^ = -I 

J^iL) ^ ah^ + PhL , (6) 

where 



h=jMO d£, and J sin^MOcos^ MOd^ ■ (7) 

Clearly, the first term in Eq. (^, proportional to 1/L, is due to the exchange interaction measured by a and the 
second term, which is proportional to L and the constant /3, represents the contribution of the magneto-crystalline 
anisotropy energy for cubic systems. Furthermore, J-{L) reaches its minimum at 




^BW-W-T7^. (8) 



Given that Ii and I2 are simple dimensionless numbers which parameterize the profile <j)o( j^), this result agrees well 
with the exact soliton solution for I1/I2 — tt^ . 

Evidently, the above discussion suggests a rather straightforward strategy for first-principles calculations: 

a. ) Assume a profile 0o(x')- 

b. ) Carry out a density functional calculation for the magnetization M{z) constrained to follow the profile 0o(x) 

and calculate the total energy or grand potential n{L). 

c. ) Minimize ^1{L) with respect to L. 

For the profile we shall always take the particularly simple choice 

r c<o 

MO = { ^7^/2 < e < 1 (9) 

2 -. 

which is compared to the soHton solution in Fig. 2. For this function /i = and I2 — ^ and hence if Q.{L) is given 
by Eq. (§ 
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which imphes ^bw — v^^'y^- Note that this Bloch Wall thickness is V2 times larger than that deduced from Eq. (||). 

Thus the Bloch Wall thickness is to some extent an ill-defined concept in the present 'prescribed profile' approach to 
the problem. As will be clear presently, our calculated ri(L) can be fitted by the above functional form very accurately 
and hence it determines the coefficients a and /? from first principles. 

For orientation we note that via the Landau-Lifshitz equation, Eq. (|l|), the free-energy functional given in Eq. (Q) 
implies a spin-wave dispersion relation. In the long wave-length limit this yields 

4/3^B , 4a/iB 2 /I IN 



thus, the spin- wave stiffness constant 



and the usual cubic (fourth-order) anisotropy constant 

^4 - y/3 . (13) 

As expected, in the presence of magneto-crystalline anisotropy the spin-wave spectrum is gaped. 

Without any account of dipolar interactions, in non-relativistic quantum mechanics the orientation of the magne- 
tization is independent from the crystal axis. As is clear from the above discussion, under this circumstance (/? = 0) 
there is no Bloch Wall of finite width. Thus, our calculation of Q{L) must include spin-orbit coupling which is the ma- 
jor source of magneto-crystalline anisotropy in transition metals p2]. Rather than treating the spin-orbit interaction 
in perturbation theory, we use a fully relativistic spin-polarized Density Functional Theory |19 20|]. 



The natural variable in the microscopic theory is the magnetization averaged over a unit cell 

= / (frrn{f) , (14) 



where Vi is the volume of the ith unit cell. We may then define the local orientation by 

ei=mjmi, (15) 

where rrii denotes the length of the vector . It is this effective local orientation that we wish to constrain to follow 
the discrete version of the profile in Eq. (^). Namely, we assume that Bp is the same within an atomic plane with 
position Zp and take 



x^p 

where for a 90° Bloch Wall, which is N layer thick 



cos 



4>pX + sin (j)py , (16) 



TT 



^p=P^ iP^h---,N). (17) 

Thus, we have to solve the Kohn-Sham-Dirac equation of DFT for the circumstance where the orientation of the 
magnetization, 0p, is changing from layer to layer within a slab of N atomic layers and is uniformally and tt/2 to 
the left and right, respectively, of the slab as shown in Fig. 2. Evidently, the relativistic spin-polarized Screened- 
KKR method |2^, which treats just such geometries, is ideal for tackling this problem. Fortunately, the code which 
implements this method scales linearly with N and, as we shall demonstrate presently, can handle 400-500 layers with 
readily available computer power. 



III. THE COMPUTATIONAL METHOD 



We performed calculations proposed in the previous section by using the spin-polarized relativistic screened 
Korringa-Kohn-Rostoker (SPR-SKKR) method [Q. Although the method is by now well-established (see also Refs. 
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Ig^,^), for completeness and since we also made developments specific to the Bloch Wall problem, we briefly outline 
our computational strategy. For an ensemble of individual scatterers the multiple scattering theory (MST) yields the 
one-electron Green function for an arbitrary complex energy z and real space coordinates f, r' 

C?(z;f,r') = ^Z5(z,f„)r5^,(z)Z5',(z,r4)t- (18) 
QQ' 



-Snnr ^ { Jq (^, r„) (z, r^J 1 9 (r„ - r'J + Z^{z,rr,)J^{zyJe{r'„ - r„)} 
Q 



where n and m label two specific sites of position vectors i?„ and R'^^, while Q and Q' denote pairs of angular 
momentum quantum-numbers {n, fJ.) and (k',/z'), respectively. The functions ZQ{z,f„) and jQ{z,rn), regular and 

irregular at r„ = r — i?„ = 0, respectively are properly normalized solutions of the Kohn-Sham-Dirac equation related 
to a single, finite-ranged potential well for which we now assume a spherically symmetric effective potential, V°^{r), 
and effective field, Bf^{r), pointing along a local (positive) z coordinate axis. In short, 



where 




with the usual Pauli matrices a and the two-dimensional unity matrix 12- The numerical solution of the above 
equation with a corresponding expression for the so-called single-site i-matrix was originally given in Refs. |^,^. A 
particular feature resulting from the approach they used is that, although the i-matrix has necessarily off-diagonal 
elements, no coupling between different £ values is present. 

We further simplify the problem by adopting the atomic sphere approximation (ASA) in which the volume of the 
sphere is taken equal to the volume of the corresponding Wigner-Seitz cell. Clearly, ASA deals with overlapping 
potentials which is, in strict sense, prohibited within the MST, however, mostly for inhomogeneous systems, this 
approach is widely used, since on one hand it gives a better description of the interstitial region than the muffin-tin 
approach, on the other hand it is conceptually much simpler than a full potential description. As implied above, our 
approach allows the orientation of the magnetization to vary from site to site as required in the case of a Bloch Wall 
where we keep the orientation to be constant within the atomic planes but let it rotate from layer to layer around 
the global z direction. Therefore, in each layer there will be a rotation R, corresponding to 9 = ^ and (j) = 4>p (see 
Eq. (p^)), which relates the f-matrix in the local to that in the global coordinate system as follows 

tgiob{^^) ^ D{R)koA^)D(R)^ , (20) 

where D.{R) denotes a matrix containing, block-wise, the irreducible representations of R. 

As the solution of the above single-site problem is relatively easy and universal in applications of MST, the evaluation 
of the scattering path operator (SPO), Tqq, (z) in Eq. (|l8|), is, essentially, the main difficulty in such calculations. 
The geometrical arrangement of the scatterers involved in the system is put, in MST, into the so-called structure 
constants, Gqq, (z) = Gqq>{z; Rn, R,n) (G"™(z) = {Gqq, (z)}), which for the relativistic case can be obtained by a 
transformation in terms of the Clebsh-Gordon coefficients from its non-relativistic counterpart (see e.g. Ref. |^^). By 
defining the corresponding matrices having both site and angular momentum indices 

t(z) - {r{z) Snm} , G(z) = {G""(z)} , t(z) = {r""(z)} , 

the SPO is given by the following matrix inversion 

x(z) = (t(z)-i-G(z))-^ . (21) 

A particular problem arises from the fact that G"™(z) is long-ranged, therefore, the inversion in Eq. (|2|) can not be 
directly performed. For a system with three-dimensional periodicity the problem can be exactly handled by making 
use of the lattice Fourier-transformation, which splits Eq. (|2^) into the corresponding fc-projections which has to be 
solved in angular momentum space only. The Bloch Wall problem we deal with exhibits, however, two-dimensional 
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(2D) periodicity in the {x, y) plane, while in the z direction the translational symmetry is broken due to the variation 
of the orientation of the magnetization. Thus employing two-dimensional lattice Fourier-transforms 



GP«(z;fc||) = ^e'"^ii«iiG(z;Cp + i?s|,4) , (22) 

where p and q denote atomic layers, generated by Cp and G^, respectively, i?|| are 2D lattice vectors and fc|| is a vector 
in the first 2D Brillouin-zone (BZ), and the new matrix notation in terms of layer indices 

t{z)^{f{z)6p,}, G(z;^||) = {GP«(z;fc||)}, t(z; fc||) - {Tf«(z; fc,,)} , 

one can write 

T(z;fc||)= (t(z)-i-G(z;fc||))"' . (23) 



Evaluating the matrix inversion in Eq. (^3|) is still demanding since the structure constants involved are long-ranged 
as far as the interlayer distances are concerned. In order to render this problem tractable the concept of 'screening' 
has been introduced in the middle of nineties |^^. Without going into details, a canonical transformation of the 
i-matrices and the structure constants, which leaves the Green function invariant, in terms of repulsive scattering 
potentials makes it possible to reduce the spatial range of the effective structure constants. As this transformation 
is independent of the real scatterers in the system, it is very useful to perform self-consistent calculations, since it 
has to be performed for the structure constants only once at the beginning of calculation. In terms of the 'screened' 
quantities Eq. ( p3| ) has exactly the same form. However, since Gj"'{z;k\i) is now well-localized, it can be truncated 
for |p — g| > n at a given n (~ 3 for FCC and BCC principal facets), which in turn implies a block-tridiagonal form 
for the matrix G(z;fc||). This blocks are related to the so-called 'principal layers' containing n subsequent atomic 
layers. By splitting our system into left and right perfect semi-infinite subsystems, in each of them the scatterers are 
all identical (bulk), and into a central region, where the potentials as well as the orientation of the magnetization can 
vary, the projection of the SPO, Eq. (p3[), onto the central region can be calculated exactly, i.e. taking into account 
all the scattering events to the left and right semi- infinite regions |^^. A remarkable feature of the method is that, if 
only the layer diagonal blocks of the SPO need to be calculated, it scales linearly with the size of the central region, 
namely TV p^ . Thus it opens the way for investigating systems with inhomogeneities extending much beyond the 
atomic scales such as Bloch Walls. 

Let us now turn to the task of performing the Brillouin-zone integration 

t"™(z)=/" d2/c|| e"'^ii^-^ii-^ii'rf9(z,fc||) , (24) 



where i?„ = Cp + i?|| and Rm = Cq + while flsz denotes the volume of the 2D BZ. In Ref. |23| we described a 
method for reducing the demand of the above BZ integration using the symmetry operations of the underlying lattice. 
This is obviously useless for the present case of the 90° Bloch Wall, since the direction of the magnetization rotates 
from layer to layer, say, from the x axis to the y axis, therefore the Bloch Wall itself is not invariant under any of the 
symmetry operations of the Civ group characteristic to the BCC(OOl) BZ. 

However, it is still useful to note that the magnetization direction described by the angle ^(z) in a 90° Bloch 
Wall perpendicular to the (001) direction of a BCC lattice and satisfying the boundary conditions </)(— oo) = and 
(/)(oo) — ^ has the symmetry property 

m - J - ^^i-^) + I ■ (25) 

Taking into account also the symmetry of the underlying lattice this implies that the Bloch Wall is invariant under a 
180° rotation around the axis (110) which is in three dimensional space represented by the matrix 




(26) 

whereby the axis of the rotation should cross the z axis at z = 0, i.e. for which 0(0) = ^. Quite clearly, = Sbw- 
Also evidently, the 2D square BZ is invariant under Sbw- 
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To make use of this symmetry, let U_ be the unitary matrix which represents Sbw in the (k, /z)-space and p' denote 
the layer onto which a particular layer p is mapped by Sbw ■ Then for the corresponding single-site t-matrices we can 
write 

f'{z)^Uf{z)U} (27) 

By introducing 

U = {C/p,}, Upq = U5,v' (28) 
or alternatively, by taking the choice of p' — —p, 

( \ 
U = ... [/ ... (29) 
\ [/ / 

it follows that 

t(z)=Ut(z)U^ (30) 

A relationship between the 2D Fourier-transformed structure constants with arguments fc|| and fc|| = S'bw^ii can 
also be established using the transformation of the real-space structure constants 

G{z- SbwR, SbwR') = UG{z; R, R') u} . (31) 
Therefore, by using Eq. ( p2| ) we can proceed as follows 

= ^ e'^fi^ii C/ G (z; Cp + - - C,) u! = 
= ^ /;i(i5„+fl;'-flf ) c/ G (z; + 4) c/t = 

= e'^ii^iic/ GP«(z; fc;,) t/te-'^n-f^il , (32) 



where we defined = S'ewGp — Gp' and i?^ = SbwR^^, both being 2D lattice vectors. Thus, similar to Eq. (|2q), 
introducing 

U(fc||) = {^p,(fc||)} , ^p,(^||) = e^'"'^5ic/ V , (33) 
Eq. ( ^ ) can be written compactly 

G(z, fc|| ) = U(fc|| ) G(z, fc[| ) U(fc|| )^ (34) 



Clearly, the matrix U(fc||) can also be used instead of U in Eq. (30), which immediately implies the following trans- 
formation property for the SPO 

x(z,fc[|) =U(fc||)tr(z,fc||)U(fc||) . (35) 

As the above relationship makes possible to calculate the SPO only in half of the BZ when performing the integration 
in Eq. (24), the computational time and memory storage request of the computer code can be reduced by a factor of 
two. Alternatively, in the actual calculations we have made use of the above symmetry of the Bloch Wall by halving 
the computational demand of the inversion in Eq. (23), while keeping all the fc|| points in the BZ. Although, this 
latter procedure is almost equivalent to what has been described above, as it facilitates the tridiagonal shape of the 
corresponding matrix [ p9[ , its use is rather limited to localized schemes whereas the former one can be regarded to be 
quite general for calculations in Bloch Wall problems. 
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For our present purposes the one-electron Green function, Eq. (|18D, can be used to calculate several quantities of 
interest such as the electronic density of states (DOS) 

n{z) -^Im J d^rTr{g{z;f,r)) {z = e + i6) , (36) 

where the energy e is real, (5 is a small imaginary part serving as parameter of Lorentzian broadening, and Tr denotes 
trace of a matrix in the four-dimensional Dirac-space, the charge density 

p{f) = --Im I dzTr{g{z-f,r)) , (37) 
71" Jc 

where C is a semi-circle contour in the upper complex semi-plane starting at the bottom of the valence band, ee, and 
ending at the Fermi energy, ep, and the spin-density 

rh{f) = --Im J dzTr{j3tg{z;r,r)^ . (38) 

Furthermore, the charge, the spin-moment and the band-energy can be obtained straightforwardly as 

dz n{z) = / d^r p{f) , (39) 



d-^r rh{f) , (40) 

and 

E^= I dsen{e) ^ --Im I dz z I d-^r Tr {g{z;f,r)) , (41) 



respectively. From Eq. ( |18D it is quite obvious that the above quantities can readily be resolved into components with 
respect to cells (layers) as we shall show them when we present are results in the forthcoming sections. 

By solving the Poisson-equation for the electrostatic potential and employing Spin-Density Functional Theory in 
the Local Approximation for the exchange-correlation potential and exchange field, self-consistent calculations 
can be performed. First we carried out self-consistent calculations for BCC bulk iron, where we have used a theoretical 
lattice constant of a=5.204 a.u. derived by recent careful full-potential calculations In that calculation we kept 
the orientation of magnetization along the (100) easy axis and used 91 fc-points in the irreducible wedge of the 2D 



BZ (IBZ) for the BZ-integrations, Eq. (24). For the energy integrals, Eqs. ( |39D and ( [40| ) we used 16 points along 
the semi-circle sampled according a Gaussian-quadrature, while the corresponding summations in angular momentum 
space were subject to a cut-off of •^max=2. We converged the Fermi energy to get the corresponding charge (Z=27) 
to an accuracy of 10~* electrons. Note, that in our present approach, because of the equilibrium with the left (and 
right) semi-infinite regions, this Fermi level has to be also used in the calculations of the (finite-size) Bloch Wall. 

As mentioned earlier, we have made use of different approximations to calculate the free-energy. The first, computa- 
tionally less demanding, one is based on the frozen-potential approximation frequently used also in magneto-crystalline 
anisotropy calculations (see Ref. and refs. therein). Briefly, we used the self-consistent bulk potential and effective 
field magnitude in each layer of the Bloch Wall and we set the orientation of the exchange field by rotating succes- 
sively the bulk t-matrix from layer to layer, Eq. (|20|), according to the prescribed function (f>o{z). In this case only one 
iteration, to calculate the single-particle (band) energy and the charges, were carried out. As the charge neutrality 
(number of particles) is not preserved for the fixed volume of the central region, the grand-canonical potential has to 
be considered. For any thickness N of the Bloch Wall, taking always the difference with respect to the ferromagnetic 
state, this is approximated by 

An{N) = AEk - ep AQ = ^ {AE^^p - ep AQp) . (42) 

p 

The other approach is a fully self-consistent one. Unfortunately, the computer power available to us was sufficient for 
only a few such calculations. Therefore, only for the case of = 60 shall we present and discuss the changes in the 
electronic structure and magnetic moments due to the Bloch Wall. 



8 



IV. RESULTS 



A. The frozen potential calculations of Bloch Wall energies 

We have calculated the Bloch Wall formation energy Ari(7V) for various values of N using the frozen potential 
approximation. Note, that the summation in Eq. has to be, in principle, taken over all the layers in the system. 
Our calculations show that layers more than about ten layers away from the wall do not contribute significantly to 
the sum although AEp and AQp by themselves differ considerably from zero even far away from the wall. 

As well-known, in the phenomenological theory the exchange energy and the magneto-crystalline energy con- 
tribute equally to the Bloch Wall energy. Experimentally, for BCC Fe the magneto-crystalline anisotropy constant, 
K4, is found to be 0.3 /iRyd per atom while careful (non-orbital-polarized) LDA calculations predict 0.1-0.2 
/iRyd (see Ref. and refs. therein and note that for a cubic system K4 — 3[i?(lll) — E{100)].) This, on one hand, 
implies, that the Bloch Wall energy normalized to one layer is expected to be of the same order of magnitude as the 
anisotropy energy and, on the other, we have to calculate the Bloch Wall energy to the same accuracy as is necessary 
in magneto-crystalline anisotropy calculations. 

The main difficulty arises from the fact that, in particular, close to the real energy axis one has to sample a high 
number of fc|| -points when performing the BZ-integration in Eq. (p^). To reduce this problem we smoothened the 
energy integrals in Eqs. (^9|) and ( ^l|) by the Fermi-function at a finite temperature T, taking into account the poles 
of it below the contour C. Both from an analysis of the integrand and by checking it numerically, it turns out that the 
contour can be deformed to infinity in the upper complex semi-plane. Moreover, only a finite number of Matsubara 
poles, Zj = eF-l-«(2j-|-l)7rA:Br (j = 0, 1, 2, . . .), has to be considered. Assuming a quadratic T dependence of An{N; T) 
due to the Sommerfeld expansion, in order to perform extrapolation to T— it was necessary to take two different T 
values only. 

In our calculations we choose 300K and 150K for these two temperatures by using 32 and 40 Matsubara poles 
with 1275 and 2926 fcy-points in the 2D IBZ for zq (3 and 1.5 mRyd), respectively. These values were shown to 
be sufficiently high to yield converged bulk anisotropy energies which we calculated to check the reliability of our 
numerical evaluations with respect to other methods and also to compare to the value that can be deduced from the 
Bloch WaU energies (see Eq. @. For T=300K and 150K we got Ki=0.U2 and 0.154 yuRyd, respectively, which 
were extrapolated to ^^4=0. 158 /iRyd at T=0. Obviously, this is a very good agreement with the results of other 
first- principles calculations. 

Let's now turn to the results for the Bloch Wall energies, associated with an area of size a^, where a is the lattice 
constant of our BCC lattice, as a function of thickness L, measured in units of a (2L/a = TV), shown in Fig. || for 
T=150K and T=300K as diamonds and crosses, respectively. The first thing to note is that the calculated points are 
very well fitted by the expression Eq. (|lO|) as drawn by solid and dashed lines, respectively. The resulting values for 

a = and p = 

4 ^8 

are listed in Table I together what fitted for T=0. 

The fitted curves have their minima at 424 and 401 lattice parameters for T=300K and T=150K, respectively. Note 
that there is a point on the T = 300K curve which is beyond the minimum and hence we can be said to have crossed 
the scale gap. The minima are very shallow due to the smallness of the anisotropy energy, thus, this is likely to be a 
feature of calculations for all cubic systems. 

From the extrapolated values of a and /3 to T = 0, A'4=0.16 /iRyd and D=262 meVA^ can be derived via Eqs. ( p^ ) 
and dl^)). By using the expression of £bw from the phenomenological theory we also obtained a thickness of 394 lattice 
parameters for the Bloch Wall. Several experimental values can be found for D, e.g. 314 meVA^ by Stringfellow 
[H, 281 meVA^ by Collins et al. ||l or 280 meVA^ by Mook and Nicklow So our value is only slightly 

below the experimental findings. The same is true for the comparison with other theoretical results found from bulk 
calculations ll5UTq] . One source of discrepancies between our and other calculations is certainly the use of different 
lattice constants. It is, however, worth noting that by using the same computer code and a spin-flip technique, in the 
scalar- relativistic limit, a value of D—300 meVA^ was calculated |3^], in better agreement with those calculated by 
others. This suggest that another source of the deviation of our present value from those of others arises from the 
relativistic approach we used, that is, the spin-orbit coupling gives not only rise to the magneto-crystalline anisotropy, 
but to some extent influences the spin-spin interaction parametrized by D. 

The value of the anisotropy constant is of greater concern. As it is generally the case in LDA calculations (l^Jl^ , the 
experimental value, 0.3 /iRyd, is roughly a factor of two bigger than the theoretical results. However, the important 
point here is that we find the same value for as given above when we determine from bulk calculations 
demonstrating the internal consistency of our calculations. 
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As implied by Eq. ( p2| ) our calculation allows to resolve the Bloch Wall formation energy into contributions related 
to layer. This than can be compared to the free-energy density of the Ginzburg-Landau theory [^. Taking the mean 
value for each layer and in terms of the parameters introduced above, for a linear 90° Bloch Wall this free-energy 
density reads 

, = f ^+/3sin2(0p) if \<p<N .43. 
^ 1^ anyway 

In Fig. |4| the solid line and crosses in the inset display the layer-resolved contributions to the Bloch Wall energy as 
calculated by our first-principles method for a 800 layer thick wall and T=150K, while the dashed curve and the 
diamonds in the inset are evaluated by using Eq. (^3|) with the the corresponding fitted parameters, a and (3 (see 
Table I.), that is no further fitting have been used. Apparently, the coincidence of the two curves is nearly perfect 
even in this 'atomic scale' resolution of the Bloch Wall energy. As is clear from Eq. (^3|), in the phenomenological 
theory the exchange contribution to the layer resolved energy is constant, which equals 2a/N'^ , inside the wall and 
zero outside. A characteristic deviation from that behavior is found for the first-principles values near the edge of 
the wall. The exchange contribution to the first layer outside the Bloch wall, i.e., to layers numbered by and 801 
in Fig. U, is exactly half of the above constant value (see inset of Fig. ^). This is due to the non-local nature of the 
exchange couplings not accounted properly within the phenomenological theory. We note that Bloch Wall formation 
can be discussed also by using an effective Heisenberg Hamiltonian with a cubic anisotropy term, in which the above 
feature is readily recovered. Small oscillations can also be seen near the edge of the Bloch Wall. These are most likely 
related to Friedel-type oscillations, which in homogeneous systems arise due to any imperfections, as relaxations in 
the electronic structure. 



B. Selfconsistent calculation 



In order to access reliably the electronic structure in the presence of a Bloch wall selfconsistent calculations have 
to be performed. These calculations are not ground state calculations because the magnetic moment in every layer 
is forced to point along a prescribed direction. The proper way to perform such calculations is to introduce a 
constraining field which forces the moment to point along the chosen direction pl| . For our first attempts we ignored 
the constraining field and, instead, took the projection of the magnetic moment onto the prescribed direction after 
every iteration. 

We have performed this calculation for a 60-layer Bloch wall with a magnetization profile accor ding to Eq. (|). 
We allowed 21 layers on the two sides outside the Bloch wall to relax. That is, given the symmetry of the Bloch 
Wall, 51 different potentials were involved in this calculation. The charging due to the Bloch Wall per atom for the 
different layers works out to be smaller than 10~^e and, therefore, can be considered to be zero within the accuracy 
of our calculations. This means we observed no charge redistribution (transfer) across the Bloch Wall. Despite of 
this fact, as inferred from Fig. ||, the densities of states show characteristic changes for different layers in the Bloch 
Wall. Interestingly, the biggest changes are found at energies with also big peaks in the bulk DOS, which presumably 
indicates lifting of some degeneracies, related to the cubic symmetry, due to the presence of Bloch Wall. 

Contrary to the charges, the magnetic moments display a small but clear deviations from their bulk value. In 
Fig. H we show the deviation of the moment from the bulk value. Obviously, the moments decrease gradually when 
approaching the center of the wall. Surprisingly, however, the moments outside the wall still quite differ from their 
bulk value, approaching it relatively far from the wall only. Similarly, at the center of the wall one would expect a 
moment equal to that calculated for a bulk with (110) hard axis. This value is shown by dashed line in Fig. 3. 

The component of the moment parallel to the wall but perpendicular to the exchange field is shown in Fig. 7[ (One 
should note that the frame of reference is a local one which turns round with the exchange field.) Not surprisingly, 
the strong peaks appear at those layers which are just outside the Bloch Wall, as these are the layers with the most 
asymmetric neighborhood. Clearly, these peaks are an artifact due to the prescribed magnetization profile which 
display pronounced kinks at the two borders of the Bloch Wall (see Fig. 2). The component perpendicular to the 
wall is smaller than 10~^/iB and therefore negligible. 

These calculations were performed with 91 fc|| -points in the two dimensional irreducible Brillouin zone. To achieve 
convergence roughly 100 iterations have been necessary (Where one iteration took roughly 14 min on 16 nodes on a 
T3E). Clearly, the number of selfconsistent iterations necessary to obtain reasonable convergence goes up with the 
number of layers and we have not been able to achieve convergence for thicker Bloch walls within the limited CPU time 
available to us. Due to the relatively small number of -points used in this calculations wc do not claim convergence 
for our results. But we expect the results to be qualitatively correct. 
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V. CONCLUSIONS 



In short, we have presented a first-principles, that is to say parameter-free and yet materials specific, description of 
Bloch Walls in ferromagnets. As an example we have prescribed the orientation of the magnetization density, M{f), 
to evolve from the easy axis (100) to (010) in BCC iron and calculated its energy, Ebw, as a function of the width, 
^BW) of the transition region. To do this we used fully relativistic spin-polarized Density Functional Theory (DFT) 
and solved the Kohn-Sham-Dirac equation by the SPR-SKKR method ||2^. Due to the relativistic description of the 
electrons, their spin and orbital degrees of freedom were treated on equal footing and, hence, the calculation gave 
a full account of the magneto-crystalline anisotropy. Consequently, the Bloch Wall energy had a minimum. The 
equilibrium width as well as the full curve, -E'bw(^bw) displayed in Fig. 3. was found to be in good agreement with 
available experimental data. 

The novel feature of this kind of electronic theory of a Bloch Wall, as compared to its conventional phcnomcnological 
description Q , is that it provides an account of the distortions of the electronic structure due to the presence of the 
Bloch Wall as well as its shape, width and energy. In particular, we have calculated the variations of the size of the 
magnetic moment, related to the local exchange splitting, and the local densities of states from layer to layer across 
the transition region. These results are displayed in Figs. 5. and 6. Evidently, these changes are small as expected 
on the grounds that BCC Fe is a good moment system. Nevertheless, even in this case they contain the essential 
information to describe the scattering of electrons by Bloch Walls as one needs to do in a study of magneto-transport 

0-0 • 

Finally, we note that our calculations scale linearly with the number of layers, N, within the Bloch Wall. Due to 
this fact we were able to perform calculations for up to 800 layers. That is to say, for the first time, we were able to 
describe a mesoscopic magnetic defect in fully first-principles terms. 
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T [K] 


a [mRyd] 


/3 [^Ryd] 


300 


64.1 


0.0356 


150 


63.2 


0.0392 





62.9 


0.0404 



TABLE I. Parameters derived from a least square fit of tlie data in Fig. 3 for T=300K and 150K to the function Eq. (^. 
In the last row the corresponding values from a quadratic interpolation to T=0 are found. 




FIG. 1. Sketch of the spherical coordinates, 9 and (p, as used throughout in the paper. 
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FIG. 2. Comparison of the soliton solution, Eq. (g), (solid line) and the linear magnetization profile, Eq. (dashed line) 
used in our calculations. The numeration of the layers and the partitioning of the system into different regions as used in the 
actual calculations is also depicted on the picture. 
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FIG. 3. Bloch wall energy as a function of the Bloch wall thickness. The diamonds and crosses stand for the results from the 
frozen potential calculations for T=150K and T=300K, respectively, with corresponding fits to the function Eq. (jl^) displayed 
in order by the dashed and solid lines. 
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FIG. 4. Layer resolved energies for a 800 layer Bloch Wall. The full line (crosses in the inset) show the result from the 
frozen potential calculation at T=150K. The inset shows a magnification of the border area. The dashed line (diamonds in 
the inset) show the layer resolved energy contributions according to Eq. ( ^3| ) of the phenomenological theory with parameters 
taken from Table I. 




-0.4 -0.3 -0.2 -0.1 0.1 0.2 
e-ep(Ryd) 



FIG. 5. The density of states (DOS) for bulk Fe BCC (thick line) together with their difference from the bulk for several 
layers in a 60 layer Bloch Wall magnified by 1000. In the legend, negative/positive numbers label layers outside/inside the 
Bloch WaU. 
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FIG. 6. The deviation of the magnetic moment from its value in the bulk, where it points in the easy direction (100), for 
each layer in a 60 layer Bloch Wall. The dashed line indicates the moment obtained in a bulk calculation with the moment 
along the (110) direction. 
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FIG. 7. The perpendicular component of the magnetic moment for each layer in a 60 layer Bloch Wall. 
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